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We report a two-detector measurement of the propagation speed of neutrinos over a baseline of 
734 km. The measurement was made with the NuMI beam at Fermilab between the near and far 
MINOS detectors. The fractional difference between the neutrino speed and the speed of light is 
determined to be (u/c — 1) = (1.0 ± 1.1) x 10“®, consistent with relativistic neutrinos. 
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PACS numbers: 14.60.Lm, 06.30.Gv 


I. INTRODUCTION 

A cornerstone of the theory of relativity is that there 
is a single limiting speed, the speed of light in a vacuum 
c, which cannot be exceeded. Observations of neutrinos 
from SN1987A[TH3] and accelerator experiments [H [S] 
have set limits on the difference between the speed of 
neutrino propagation and that of light, all consistent with 
V = c. In September 2011, the OPERA experiment re¬ 
ported a measurement [5] , in striking conflict with both 
theory and experiment, which has since been revised to 
resolve the inconsistency [7] . The initial OPERA news 
motivated a number of further measurements [HHE] . We 
report here a new precision measurement of the speed of 
neutrinos using the NuMI muon neutrino beam at the 
Fermi National Accelerator Laboratory (Fermilab) [T3] 
and the two MINOS detectors m, using a significantly 
upgraded time synchronization system and an exposure 
of 0.8 X10^° protons on target in March and April of 2012. 
This measurement has similar precision to the CERN to 
LNGS measurements referenced above, but with two sig¬ 
nificant differences in addition to being made in a dif¬ 
ferent lab with different dominant uncertainties. First, 
the NuMI beam has a mean neutrino energy of 2.8 GeV, 
almost an order of magnitude lower than the 17 GeV 
GNGS beam. Second, the measurement presented here 
measures the neutrino time-of-flight using two neutrino 
detectors, rather than using proton bunches in the ac¬ 
celerator to obtain the start time and neutrinos for the 
stop time. This neutrino time-of-flight measurement is 
the most precise ever, although the velocity precision is 


~9.6 /is spill, 6 batches 




18.83 ns 100 ns 0.8 ns RMS 


FIG. 1: Schematic of the spill structure showing the 6 batches 
(top), the gap between two of the batches (bottom), and the 
bunch gap and width. 
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more limited by distance uncertainties. 

The neutrino velocity measurement is conceptually 
straightforward, consisting of a measurement of the dis¬ 
tance between the two detectors and the time it takes for 
a neutrino to pass between them. We can never observe 
the same neutrino in both detectors since the process of 
detection is destructive. We address this issue by making 
two separate measurements with respect to the same time 
reference. Specifically, we first make a measurement of 
the time of arrival of a bunch of neutrinos in the MINOS 
Near Detector (ND) referenced to a time marker derived 
from the proton beam current near the neutrino produc¬ 
tion target and then of that bunch’s arrival at the MINOS 
Far Detector (ED), also with respect to the proton time 
reference. Subtraction of these two measurements, cor¬ 
rected for various offsets from the detection process, gives 
the time-of-flight of the neutrinos over the distance be¬ 
tween the two detectors. Care is taken so that when the 
subtraction of the two times is made, the major part of 
the uncertainties in the detection cancels, leaving a high 
precision determination of the time the neutrinos took to 
travel from the near detector to the far detector. 

The measurement uses a system of time transfer be¬ 
tween the beam current measurement and each detector 
with sets of periodically calibrated Global Positioning 
System (GPS) receivers and atomic clocks. Two Way 
Satellite Time and Transfer (TWSTT) is also used as 
an independent technique to calibrate the time offsets. 
The combination of these techniques allowed reliable es¬ 
timates of the time synchronization errors, yielding a very 
robust measurement that has the smallest error for the 
flight time of neutrinos ever achieved. 


II. BEAM AND DETECTORS 


The neutrino beam |I8j is produced at Fermilab by 
120GeV/c protons striking a graphite target. The re¬ 
sulting positively charged pions and kaons are focused 
by pulsed magnetic horns and then allowed to decay 
in a 675 m long helium-filled volume, producing a 
dominated neutrino beam with a peak in the event en¬ 
ergy spectrum at about 3 GeV and a tail at higher ener¬ 
gies m The time structure of the beam is illustrated in 
Fig.i During acceleration, the protons are grouped into 
six batches, each approximately 1.6 /iS long and separated 
by about 100 ns. Each batch consists of 81 bunches which 
are 0.8ns RMS wide (3.5ns full width at base), spaced 
at 18.83 ns intervals, resulting from the Main Injector’s 
53.103480 MHz synchrotron acceleration. 
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A. The MINOS Detectors 

The two MINOS detectors [ll] are steel and scintilla¬ 
tor tracking calorimeters with toroidal magnetic helds 
averaging 1.3 T in the steel. Each detector consists of 
2.54 cm thick steel plates interleaved with 1 cm thick plas¬ 
tic scintillator planes. The scintillator planes are com¬ 
posed of 4.1cm wide strips. Scintillation light is read 
out by multi-anode photomultiplier tubes (PMTs) via 
wavelength-shifting fibers. The 0.98 kton ND is located 
1.04 km downstream of the production target and 104 m 
underground. The 5.4 kton FD is approximately 735 km 
downstream of the target and 705 m underground. 

Muon neutrinos are identified in the detectors through 
their charged-current interactions, + A ^~ + X 
where A is the target nucleus (normally Fe) and X rep¬ 
resents the final state which may contain pions, other 
hadrons, and nuclear fragments. The muon typically 
leaves a well-defined energy deposit in the detector, cross¬ 
ing tens of scintillator planes that can be reconstructed 
as a muon track. Due to the high rate of neutrino inter¬ 
actions within the near detector there are differences in 
the way the times of the individual pulses are recorded 
in the two detectors. The subsequent reconstruction and 
selection of the interactions, however, are done in an iden¬ 
tical way. The time and position of the neutrino interac¬ 
tion vertex is calculated using the TDC (time-to-digital 
converter) times of the PMT pulses, combined with the 
known spatial geometry of the scintillator strips produc¬ 
ing the light. 

The resolution on the time of the neutrino interaction 
is 1.5 ns within both detectors. To achieve this, the detec¬ 
tors are calibrated by first applying a “time-walk” time- 
slew correction for the variation of time with pulse height 
and then a strip-by-strip offset obtained by minimizing 
the residual offset in an ensemble of muon tracks. The 
FD time resolution capabilities and calibration |16j were 
established to distinguish upward atmospheric neutrino 
events from downward cosmic rays verified in the 
CalDet test-beam experiment CHI, and the detailed cal¬ 
ibration constants were obtained with cosmic ray muon 
events collected over a number of years. The ND cali¬ 
bration was performed in a similar way using secondary 
muons from the beam and was cross-checked with cosmic 
rays. The time resolution of both detectors is sufficient 
to resolve the bunch structure shown in Fig. and was 
cross-checked with portable counters [TS]. A study of the 
broadening of the spill structure by the detector resolu¬ 
tion was made at the ND by looking at the reconstructed 
times of neutrino interactions within the same accelerator 
spill to resolve the beam time structure. 

The MINOS experiment has previously reported a 
measurement of neutrino speed in Ref. [S]. Since then, 
the experiment has collected a factor of 8.5 times more 
data. Additionally, a comprehensive study of the com¬ 
ponents of the original MINOS timing systems was con¬ 
ducted and several new corrections have been applied to 
this larger dataset. Some of these studies required com¬ 


ponents of the new timing system, described below, to 
be operated in parallel with the original system; or by 
specific tests, comparing measurements between signals 
in the new and original systems. In particular, it was 
found that a random offset on the order of 20 ns was intro¬ 
duced each time the original GPS receivers were powered 
on, which is not unknown in GPS receivers not specifi¬ 
cally designed for precise time. An offset of this size 
was covered by the systematic errors quoted in [5] . The 
analysis of the full MINOS data sample, using the origi¬ 
nal MINOS timing system and new corrections, yields a 
systematic error dominated fractional neutrino speed of 
(v/c-l) = (0.6± 1.3) X 10-^ 

The order-of-magnitude more precise measurement de¬ 
scribed in this paper uses a new timing apparatus, mak¬ 
ing it insensitive to the variations of the old system. The 
rest of this paper describes this new system and a new 
analysis based only on the data taken after the new sys¬ 
tem was installed. 


B. Proton Beam Measurement 

The time profile of the proton beam is measured using 
a resistive wall current monitor (RWCM) [5D] situated 
along the beam pipe, which is between the extraction 
point from the Main Injector and the NuMI target. 
The RWCM consists of a resistive network bridged across 
an electrically insulating ceramic break in the stainless- 
steel pipe. As the beam passes, an image current is in¬ 
duced on the surface of the pipe, creating a measurable 
voltage across the resistive network. The voltage signal 
from the device is measured for each spill with a wave¬ 
form digitizer with 1.5 GHz analog bandwidth, which is 
the limiting bandwidth of the system. 


C. Event Timing 

The data presented here use timing components shown 
in Fig. Original local timing systems used to internally 
synchronize different parts of each detector are retained, 
as they are integral to the experiment’s data acquisition 
system. The new timing system is used to time stamp the 
old system’s timing markers, allowing more precise calcu¬ 
lations of neutrino interaction times without disturbing 
the well-tested and robust means of acquiring the data. 
The new system’s GPS units and its overall synchroniza¬ 
tion are significantly upgraded, as shown in the upper 
part of Fig. Stable atomic reference clocks are in¬ 
stalled at each detector location. The manner in which 
timing synchronization is transferred between the surface 
and the underground detector locations is also upgraded 
with optical fibers operating in both directions, transfer¬ 
ring both I Hz (or pulse-per-second, PPS) and 10 MHz 
signals and allowing continuous monitoring of the delays 
in the links. 
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FIG. 2: Simplified layout of the main time synchronization components of the experiment. The beam and detectors are shown 
across the middle of the diagram, the newly added time measurement apparatus is shown above the beam and original time 
synchronization (which is still used) is shown underneath the beam. Points marked “T” show where an interval timer is 
used permanently and points marked “X” show a connection point where equipment is connected for short periods during 
the run (the traveling GPS units and the TWSTT (Two-Way Satellite Time Transfer) equipment). Many of the links are 
synchronized/monitored with both PPS and 10MHz signals using TWTT (Two-Way Time Transfer). Not shown are the 
portable detectors. 


A local Cs atomic clock is installed at each detector, 
and a Rb clock is installed at the RWCM site to estab¬ 
lish free running PPS and 10 MHz time references at each 
site. Distribution amplifiers (high output-to-output iso¬ 
lation and low-jitter signal fan-outs) are used to define 
the reference points and provide multiple ports to make 
timing measurements. The internal time synchronization 
of each detector is measured with respect to these tim¬ 
ing reference points using interval timers. A number of 
these interval timers are installed permanently in the ex¬ 
periment and used for the duration of the time-of-flight 
experiment as shown in Fig. The offsets between the 
three atomic clocks at the different sites are continuously 
measured with the new system of GPS receivers and in¬ 
terval timers. 

At the FD, the time synchronization within the de¬ 
tector uses a 40 MHz clock and PPS boundary markers 
derived from the original FD GPS receiver, which runs 
independently of the new timing equipment. The PPS 
are encoded on the 40 MHz clock signals to allow distri¬ 


bution to the 46 readout boards in the detector over a 
single network of same-length cables. The times of the 
scintillator pulses in the detectors are measured using 
a TDC implemented by multiplying the 40 MHz to give 
160 MHz, then using four delayed versions of this clock 
to generate an effective 640 MHz clock. 

The offset between the new local Gs clock and the PPS 
from the original FD GPS receiver (used for the internal 
detector time synchronization) is measured and recorded 
each second with an interval timer. For each neutrino 
interaction, for which data acquisistion is started using a 
simple activity trigger, we correct the measured detector 
time to the new local Cs clock time using the interval 
timer measurement closest in time to the neutrino event. 
Variations in the frequency of the 40 MHz FD clock are 
monitored by measuring the interval between successive 
PPS signals with the Cs clock which is observed to vary 
smoothly by typically ±10 ns in the short term. As a 
cross-check of the stability of the timing distribution sys¬ 
tem, signals were read from one of the readout boards 
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and compared with the new time reference system using 
an interval timer. 

At the ND, the acquisition of data is started by a trig¬ 
ger signal from the accelerator which arrives about 20 fis 
before the neutrino beam. The arrival time of this trig¬ 
ger signal within the ND time distribution electronics is 
measured against the local Cs clock reference point with 
an interval timer. The ND measures the time of the PMT 
pulses from neutrino interactions relative to the trigger 
signal using a local 53.1 MHz clock from a crystal oscilla¬ 
tor that is distributed to the readout boards at the ND. 
Cross-checks were also made with interval timer measure¬ 
ments between points in the old and new timing systems. 

The 53.1MHz clock was found to have a 100 Hz vari¬ 
ation synchronous to the accelerator cycle. However, 
this clock is only used to measure short time intervals 
of O{20fj,s), so the effect is around 40 ps at worst and is 
thus neglected. 

Similarly, at the RWCM, the same 20 /iS accelerator 
signal used to start digitization at the ND is also used to 
trigger the RWCM digitizer and is measured against the 
Rb clock with an interval timer. The early accelerator 
trigger signals at both the ND and RWCM are derived 
from the same source, although this is not necessary for 
the time-of-flight measurement. This is simply a signal 
which precedes the arrival of the neutrinos by a small 
amount, used to trigger the detectors, and measured with 
respect to the local timing reference point by the interval 
counter. 


D. Detector Time Offset 

The mechanisms by which the detectors operate cause 
systematic time delays (latencies) between the neutrino 
interaction time in the detector and the time when it is 
recorded. A separate portable detector is used to mea¬ 
sure the latency of each MINOS detector. The portable 
detector consists of a pair of planes, each of active area 
63 cm X 57 cm constructed from eight plastic-scintillator 
strips left over from the original MINOS construction 
plus eight similar 3 cm wide newly constructed strips. 
The two parallel planes are stacked and oriented so their 
strips are rotated by 90° from each other. Coincident 
signals from the two planes are timed with respect to the 
PPS from the local Cs clock reference point. 

The portable detector was first placed immediately be¬ 
hind the ND where some of the muons created in the ND 
volume by beam neutrinos could pass through it. By 
matching PMT pulses in time and correcting for longi¬ 
tudinal position, we obtain a relative latency measure¬ 
ment between the portable detector and the ND. This is 
measured to be (36 ± 4) ns. Following this, the portable 
detector was transported to the FD and by using cos¬ 
mic ray muons which passed through both the FD and 
the portable detector, the relative latency between the 
portable detector and FD was measured to be (12±4) ns. 
Most of that up to 4 ns unknown latency comes from un¬ 


certainties on the propagation time of the signals in the 
portable detector itself: but the same detector and elec¬ 
tronics were used at both ND and FD. Thus, when com¬ 
paring the ND and FD latencies quoted above, the un¬ 
known delays in the portable detector itself subtract out, 
leaving a relative ND-FD latency of (24±1) ns. The error 
assigned covers the jitter on the latency measurements, 
small drifts observed over time, differences in magnetic 
fields at the PMTs at the two locations, and the some¬ 
what different energies of the different samples of muons 
used in the measurement. A second identical portable 
detector was used together with the first one at the ND 
for studies to characterize the resolution and stability of 
the counters and electronics. 


E. Detector and Baseline Survey 

The straight-line distance between the front faces of 
the near and far detectors has been determined to 70 cm 
precision. This section describes how this precision is 
achieved. 

Survey control networks have been established on the 
surface and underground at the ND and FD sites. Both 
detectors have been located relative to their respective 
underground survey control networks to within 0.5 cm. 
The two surface control networks are connected via high 
precision GPS measurements to an overall accuracy of 
about 1 cm. The tie between the surface network and the 
ND underground control network is straightforward and 
has been accomplished utilizing standard optical survey 
methods with millimeter accuracy. These measurements 
are described in Refs. [HIllllS]. 

For the FD however, there is no direct plumb line down 
the sloped shaft and issues with atmospheric stratifica¬ 
tion prohibit optical surveys. Therefore, a Honeywell In¬ 
ertial Navigation Unit (INS) containing three gyroscopes 
and three accelerometers was utilized to connect the sur¬ 
face and underground control networks. The INS was 
mounted in the elevator cage and traveled multiple times 
up and down the mine shaft, stopping each time at four 
approximately equal distance positions to reset accumu¬ 
lated velocity errors in the INS. Limiting factors of the 
accuracy include: the relatively high vibration rate of the 
elevator cage; the fact that the cage stops at slightly dif¬ 
ferent places each time; and residual oscillations as the 
cage came to a stop. The INS measurement is detailed 
in Ref. [23]. 

Observations from the gyroscopes and accelerometers 
of the INS were used to connect the FD’s underground 
coordinates to the surface and establish the detector’s 
NAD83 |25| coordinates. NAD83 is the horizontal control 
datum for North America, and the GRS80 [2S] reference 
ellipsoid was used in these conversions. The position co¬ 
ordinates of the centers of the front faces of each detector 
were both converted to these geocentric coordinates in or¬ 
der to compute the Euclidean distance between the two 
detectors. Two independent coordinate transforms were 
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done, agreeing with each other to 0.2 ns. The resulting 
longitudinal uncertainty from the front face of the ND 
to the front face of the FD of 70 cm is dominated by 
the limited INS repeatability measurements in the mine 
shaft. 

The resulting time-of-flight uncertainty caused by the 
positional uncertainties in this distance is dominated by 
the INS error and totals 2.3 ns: this is the dominant sys¬ 
tematic error in the final neutrino speed calculation. If 
we instead take the speed of highly relativistic neutrinos 
to be given as c, we can turn this measurement around to 
make a neutrino-based survey of the location of the FD. 
This interpretation of the time-of-flight data presented 
in the conclusion of this paper below (Sec. IVI suggests 
that the FD is (0.72±0.03(stat)±0.39(s?/st)) rn closer to 
Fermilab than the inertial survey indicates. 


III. TIME SYNCHRONIZATION 

The time synchronization between sites was imple¬ 
mented using several independent techniques (GPS, two- 
way fiber-based and two-way satellite-based) and also 
two different processing methods (code-based common 
view and carrier-phase based common-view for GPS 
data) [23 • This redundancy ensured robustness and al¬ 
lowed for assessment of systematics. The timing systems 
used in the experiment are shown schematically in Fig.[^ 
Two GPS receivers were deployed at each of the three 
locations (RWGM, ND, and FD) and a two-way time 
transfer (TWTT) system using fiber links was installed 
between the RWGM and ND. A two-way satellite time 
transfer (TWSTT) synchronization was also performed 
over a 36 hour period during the data taking period us¬ 
ing a dedicated satellite link between the ND and FD. 
The ND and FD atomic clock time references are un¬ 
derground and were transferred to secondary references 
on the surface using TWTT fiber links. The three fiber 
links in the experiment are similar and transfer both PPS 
and 10 MHz signals on independent sets of fibers. All the 
timing instruments (distribution amplifiers, optical link 
transmitters and receivers, interval timers) were commer¬ 
cial units specifically intended for transferring or measur¬ 
ing time accurately. Gables were high quality RG/316- 
DS using SMA or BNC connectors, and obtained specif¬ 
ically for this measurement. 


A. Fiber links 

The TWTT technique is used in the three fiber links 
shown on Fig. [^by employing a second fiber in the same 
bundle, so the fiber is therefore susceptible to the same 
environmental changes. The timing signal is sent out 
from the reference point at the first location to the second 
location where it is used to establish a second reference 
point. It is then sent back on the second fiber to the first 
location. At the first location, the delay between the ar¬ 


riving signal and the first reference point is continuously 
recorded using an interval timer and used to correct the 
data. 

The largest variability in the round-trip time is in 
the FD surface-underground link, which shows a ther¬ 
mal day/night effect of 200 ps round-trip and a similar 
overall variability during the run. This is corrected to 
better than 50 ps by continuous monitoring of the round 
trip time of the PPS along this link. The differences in 
delays of the communication fibers, transmitters and re¬ 
ceivers are determined by swapping modules. 

A portable Cs clock was used to verify the calibration 
of the surface to underground links by measuring the 
offset between portable and reference clocks, first on the 
surface, then underground, then on the surface again and 
correcting for the relative clock drift. The measurements 
show an average discrepancy of 400 ps at the ND. At 
the FD, the internal delays in the optical receivers were 
not measured directly, but corrected using the average 
of a series of three portable clock measurements made 
on different days, and the ±550 ps maximum deviations 
of the three measurements is applied as the systematic 
uncertainty of this correction. 


B. Global Positioning System 

The GPS timing infrastructure consists of eight simi¬ 
lar dual-frequency GPS receivers: six identical receivers 
(Novatel OEMV) and two newer versions from the same 
manufacturer (Novatel OEM6). All receivers use anten¬ 
nas from Novatel with Andrew FSJ1-50A antenna cables 
with small (-0.028 to ±0.036 ps/m/°G) temperature coe- 
ficients. The antenna cables were annealed for tempera¬ 
ture stability of the propagation delay before installation. 
The group delays of the antenna cables were measured at 
both the LI (1575.42 MHz) and L2 (1227.60 MHz) GPS 
signal frequencies prior to installation and confirmed us¬ 
ing time-domain reflectometry. The measured cable de¬ 
lays are used to calculate the time difference between the 
two locations. 

Two receivers are located at each of the three sites. 
The antennas are located so as to minimize “multipath” 
errors from reflected GPS satellite signals arriving at 
slightly later times. The two remaining receivers were 
transported (with their antenna and antenna cable), be¬ 
tween the three MINOS sites and the National Institute 
of Standards (NIST) in Boulder, GO, to provide multiple 
differential calibrations of the fixed GPS systems. The 
two mobile receivers were used at all the sites in differ¬ 
ent orders. They were moved at approximately weekly 
intervals, spending about three active days at each site. 
When two co-located GPS receivers are calibrated, the 
difference between the local time as measured by each 
receiver is corrected for the internal delays of the two 
receivers. The GPS timing infrastructure is described 
inHHIMl. 

At each location, the timing reference signal from the 
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MJD [days] 


FIG. 3: Comparison between TWSTT and the GPS system 
between ND and FD, during the USNO TWSTT test (April 
18-19, 2012), expressed as Modified Julian Dates (MJD) in the 
figure. In the legend, “TWdirect” refers to calibrated TW¬ 
STT between antennas at ND and FD, while “TW” refers 
to the double-difference obtained from subtracting observa¬ 
tions of ND-USNO from FD-USNO. The subscripts below 
GPS identify the individual GPS receivers used. The mean 
differences used in the analysis are shown. 


atomic clock is input to the GPS receivers. Each receiver 
makes measurements of the time offsets between the tim¬ 
ing reference signal and the time obtained from the GPS 
satellites which are visible. The GPS data were processed 
in Gommon-View (CV) mode using GPS LI G/A code¬ 
based data to compute the calibrated differences between 
the atomic clocks at the various locations [5S] . The pro¬ 
cessing of carrier phase-based data was done via PPP 
(precise point positioning) algorithms that also provide 
an accurate position of each GPS antenna through the 
corrections computed by the International GNSS Service 
(IGS). The GPS data reductions were done by the au¬ 
thors at the National Institute of Standards and Tech¬ 
nology (NIST) using Natural Resources Canada’s CSRS- 
PPP online service [30] and confirmed by the authors at 
the US Naval Observatory (USNO) using both NRCan 
and the JPL automated data reduction service, which 
has independent software (GIPSY/OASIS [3T]). These 
reductions provide a large standard set of corrections in¬ 
cluding relativistic effects and the “IGS Final” correc¬ 
tions for the orbit and atomic-clock time of the satellite. 

The CV code-based processing of the difference be¬ 
tween two co-located receivers gives an estimate of the 
hardware stability (GPS receiver, antenna, and cable). 
The time differences over the time-of-flight measurement 
period showed better than 200 ps (Icr) stability at the 
RWCM and ND. The GPS receivers and other timing 
equipment at the FD are in an environmental chamber 
that holds the temperature to ±1K. This results in a 
time stability of better than 90 ps, estimated using the 
RMS and the total time deviation [35] . 


C. Synchronization Between Sites 

An estimate of the stability and accuracy for the actual 
synchronization between sites is obtained by differencing 
data from independent time-transfer links between each 
pair of sites, and where possible between independent 
units at the same site. Effectively this is calculating a 
double-difference, since each link is already a difference 
between the pair of clocks at each location. Each tech¬ 
nique and each deployed system can have a bias due to 
a variety of causes, such as multipath reflected signals, 
temperature, humidity, and mismodeling of the atmo¬ 
sphere, satellite orbit, or ionosphere. These errors can 
occur during the initial link calibration, and can also vary 
over time. Double-differences wherein the two links are 
based on different techniques (i.e. GPS GV and TWSTT 
or GPS CV and fiber-based TWTT) are particularly use¬ 
ful in constraining the biases, since they are expected to 
be independent. 

A dedicated two-way satellite link was used by US 
Naval Observatory (USNO) personnel to measure the 
time difference between ND and FD references using 
TWSTT. This involves exchanging signals between lo¬ 
cations on a bi-directional satellite link, so that environ¬ 
mental and atmospheric effects are almost the same in 
each direction. The time difference calculated by double- 
differencing time-transfer between ND and FD with one 
GPS-based link and with the TWSTT link shows sta¬ 
bility better than 800 ps and a mean difference (accu¬ 
racy) of -480 ps. A secondary TWSTT mode was also 
employed using the USNO’s facility in Washington DC 
as a transfer point, and measuring FD-ND by the double¬ 
difference between USNO-ND and USNO-FD. No truly 
relativistic corrections are required in the TWSTT anal¬ 
ysis apart from the formalism to correct for the rotation 
of the Earth (the Sagnac correction) [32]. Figure]^ shows 
the comparison between these time transfer techniques. 

The same double difference computed between one 
GPS-based link and the TWSTT link can be performed 
between two GPS-based links, one using a pair of GPS 
receivers (one at each end of the link) and the other using 
the other pair. The comparison of the two GPS receivers 
at each location also allows the detection of eventual dis¬ 
crepancies between them. Since the calibration method 
simultaneously sets all GPS links between the two sites, 
the GPS-only double-differences constrain the calibration 
variation between calibrations. On the ND to FD link, 
the stability over the entire neutrino data collection pe¬ 
riod is better than 200ps (Icr). The mean difference of 
205 ps, suggests a systematic uncertainty of 150 ps (half 
the error attributed to each link). We apply an overall 
systematic error of 500 ps due to the uncertainty in the 
overall FD-ND synchronization during the neutrino data 
collection period. 

By using the two traveling GPS units it is possible to 
perform repeated calibrations of the stationary receivers 
at the various locations. Use of more than one calibra¬ 
tion and the presence of a second receiver at each loca- 









tion allows the determination of possible time steps in 
any of the receiver’s internal time base. As an example, 
the double-difference between ND and FD, when calcu¬ 
lated using the result of the same calibration over ap¬ 
proximately eight months, showed a step of 2 ns. When 
the values of subsequent calibrations were considered, it 
became clear that the internal time base of one of the 
stationary receivers at the ND had a step of about the 
same magnitude. Two different calibration values were 
thus used for this receiver to remove the step from the 
double-difference, eliminating it from the error analysis 
for the synchronization. 


IV. DATA ANALYSIS 

We select two categories of events whose selection cri¬ 
teria are described in Ref. [33]: contained charged- 
current (CC) events, which originate in the fiducial re¬ 
gion of each detector, and, for the FD only, partially 
reconstructed events. Partially reconstructed events are 
either CC neutrino interactions in the rock surround¬ 
ing the FD producing an entering muon, or events where 
the neutrino interaction point occurred outside the de¬ 
tectors’ fiducial region. Partially reconstructed events 
usually arrive later than contained events, as the neutrino 
produces a muon at an angle, resulting in a trigonometric 
increase to the muon path-length (the muons are highly 
relativistic and their slightly slower-than-light speed does 
not contribute significantly to the later arrival). Both 
samples are subject to a further track-fit quality cut to 
ensure that we only use events with a well-measured in¬ 
teraction time. During March and April 2012, an expo¬ 
sure of 0.8 X 10^° protons on target was collected, yield¬ 
ing 195 fully contained and 177 partially reconstructed 
events which are selected at the FD. 

The data analysis proceeds by using the measured 
time distribution of the protons in the RWCM to form 
a likelihood distribution for the time-of-flight of each 
neutrino event to each detector. A likelihood distribu¬ 
tion is required because, for a given neutrino, we do not 
know which part of the accelerator spill (as shown in 
Fig. 0 is responsible for producing the observed neu¬ 
trino. The RWCM timing waveform is convolved with 
a 1.5 ns (RMS) wide Gaussian distribution representing 
the detector timing resolution and shifted by the pre¬ 
dicted time-of-flight to form a probability distribution 
function for the contained events. For partially recon¬ 
structed events, the probability density function is fur¬ 
ther convolved with a delay distribution to take account 
of the increased muon path length calculated using Monte 
Carlo simulations. 

We multiply the event likelihoods of all the observed 
events together to obtain an overall probability as a func¬ 
tion of the time-of-flight parameter. The time-of-flight 
which gives the maximum combined probability is then 
our measurement of the overall neutrino flight time from 
the RWCM to the detector. This procedure is followed 
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FIG. 4: Daily variation in time-of-flight between the RWCM 
and ND. The GPS (red points) and TWTT (black points) 
techniques were used for transferring the time between the 
two sites. The difference between these two techniques is 
confirmation of the stability of the time synchronization. 
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TABLE I: Dominant systematic uncertainties (Tic) 


Systematic Uncertainty Value 

Inertial survey of the FD location 2.3 ns 

Relative ND-FD latency 1.0 ns 

FD TWTT between surface and underground 0.6 ns 
GPS time-transfer accuracy 0.5 ns 


at each of the two detectors. 

To obtain the final result for the ND to FD time-of- 
flight, the two times obtained are subtracted, eliminating 
any systematic offset associated with the beam measure¬ 
ment and the portable detector measurement. Having 
two neutrino detectors, near and far, is unique to this 
measurement, as other neutrino time-of-flight analyses 
must rely on beam measurements. 

At the ND, the neutrino interaction rate is sufficient 
to measure the time-of-flight for each day with statistical 
error below 50 ps, permitting a test of the long-term sta¬ 
bility of both the GPS and the TWTT timing systems 
using neutrinos directly, over the short baseline between 
RWCM and ND. The results from both these methods 
are shown in Fig. The surveyed distance between the 
RWCM and the ND is combined with the absolute la¬ 
tency of the ND measured with the portable detector to 
create an expected time-of-flight of 4622.7±4.0ns, where 
the uncertainty comes from the absolute latency of the 
portable detector and the RWCM-ND synchronization. 
Also included is a small effect of 0(500 ps) from the pro¬ 
tons and pions traveling slightly below the speed of light 
and away from the beam axis. 

Figurej^shows that the daily measurements from both 
methods are contained within a 1 ns range consistent 
with this expectation. Note that the precision of this 
test is not as good as on the final time-of-flight result 
which benefits from the cancellation of uncertainty in 
the subtraction technique. Figure]^ shows that the time- 
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FIG. 5: Arrival time distribution at the near detector modulo 
the 18.83ns bunch separation (blue histogram), which is well 
fit by a gaussian with a 1.6ns sigma (red line). The high 
neutrino statistics available at the ND establish the shape 
and timing of the neutrino bunch structure. 

of-flight measurements are stable to about 200 ps when 
using the GPS for time-transfer and to better than 50 ps 
when using the TWTT for time transfer. As the GPS 
time was common to both detectors and crosschecked by 
the USNO TWSTT measurement, the mean GPS-based 
value of 4621.1ns is used later in calculating the final 
time-of-flight, with the difference being covered by the 
systematic error described in Sec |IIIG[ 

Figure shows the neutrino arrival time distribution 
at the ND after the end of the preceding 18.83 ns long 
bunch, and Fig.|^the same at the FD. The high neutrino 
statistics available at the ND illustrate that the neutrino 
production in a bunch is well fit by a Gaussian with a 
1.6 ns sigma: this width is driven by the per-event 1.5 ns 
detector time resolution rather than the bunch width it¬ 
self. The contained event sample at the FD is consistent 
with that distribution, confirming the resolution and sta¬ 
bility of the time measurement system. The partially re¬ 
constructed events at the FD are also shown in Fig. 
and are seen to arrive later and with a bigger spread as 
expected. 

Gombining the contained and partially reconstructed 
samples, the time-of-flight between the RWGM and FD 
is found to be (2 453 935.0 ± 0.1) ns, considering only 
statistical errors. Subtracting the measured time-of- 
flight (using GPS) between RWGM and ND of 4621.1ns 
we obtain the time-of-flight between ND and FD as 
(2 449 313.9 ± 0.1) ns (statistical error only): the most 
precise measurement of the neutrino time-of-flight ever 


achieved, and the only one obtained directly using two 
neutrino detectors. The time required to traverse the 
distance between the front face of the Near and Far de¬ 
tectors at the speed of light, including the Sagnac cor¬ 
rection, is (2 449 316.3 ± 2.3) ns, where the dominant un¬ 
certainty comes from the inertial survey of the FD loca¬ 
tion. Gombining these, together with the other sources 



FIG. 6: Arrival time distribution at the far detector mod¬ 
ulo the 18.83 ns bunch separation for contained (blue hatched 
histogram) and partially contained (red line) events. The ap¬ 
pearance of the bunch shape in the FD demonstrates that the 
timing system is functioning as expected. 


of systematic error listed in Table |Ij yields a value for the 
difference in arrival time of the neutrino and the speed 
of light prediction of d = (2.4±0.1(stat.) ±2.6(s?/st.)) ns. 
The fractional neutrino speed is therefore found to be 
(yjc — 1) = (1.0 ± 1.1) X 10“®, consistent to 1 cr with 
relativistic neutrinos. 
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